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Abstract 

We report preliminary results on time-dependent CP asymmetries in B —> D*^^ decays. The 
CP asymmetry in these decays is proportional to 2Rd*-k sin(20i + 03 ± 5d*k), where Rd*-k is 
the ratio of the magnitudes of the doubly-Cabibbo-suppressed and Cabibbo-favoured amplitudes, 
5d*tt is the strong phase difference between them, and 0i and 03 are two angles of the CKM 
Unitarity Triangle. This study is based on a large data sample that contains 657 million BB 
pairs collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider at the T(45) 
resonance. We use a partial reconstruction technique, wherein signal B — > L'* =F 7r =l= events are 
identified using information only from the tt^ from the B decay and the charged slow pion from 
the subsequent decay of the D*~ . We obtain S + (D*n) = +0.057 ± 0.019(stat) ± 0.012(sys) and 
S-(D*ir) = +0.038 ± 0.020(stat) ± O.OlO(sys). 

PACS numbers: 
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INTRODUCTION 



In the Standard Model (SM), quark flavour mixing occurs via the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [1]. CP violation in the SM occurs due to the presence of a complex 
phase in the CKM matrix. Precision measurements of the parameters of CKM matrix are 
of utmost importance to constrain the SM and measure the amount of CP violation. The 
study of the time-dependent decay rates of B°(B°) — > D* T 7r ± provide a theoretically clean 
method for extracting sin(20! + 3 ) [2], where 0i and 03 are angles of the CKM Unitarity 
Triangle. As shown in Fig. 1, this decay can be mediated by both Cabibbo-favoured (CFD) 
and doubly-Cabibbo-suppressed (DCSD) processes, whose amplitudes are proportional to 
VdXud and V* b Vcd respectively, which have a relative weak phase 03. 
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FIG. 1: Diagrams for B° -> D*~tt+ (left) and B° -> D*~tt + (right). Those for B° -» D*+tt~ and 
B° D*+ TT can be obtained by charge conjugation. 

The time- dependent decay rates are given by [3] 

P(B° -> £) w± 7r T ) = —e~ lAtl/T B x [l =f C cos(AmAt) - S' ± sin(AmAt 

P(B° -> j D ( * )± tt t ) = — e - lAtl/T z° x [1 ±Ccos(AmAt) + S ± sm(AmAt 

8t b o 

Here At is the difference between the time of the decay and the time that the flavour of the 
B meson is tagged, r^o is the average neutral B meson lifetime, Am is the B°-B° mixing 
parameter, and C — (1 — B 2 ) / (1 + R 2 ), where R is the ratio of the magnitudes between 
the DCSD and CFD (we assume the magnitudes of both the CFD and DCSD amplitudes 
are the same for B° and B° decays). The CP violation parameters are given by 

± _ 2(-l) L fisin(20 1 + 3 ±^) 

b ~ oTrJ) ' { ) 

where L is the orbital angular momentum of the final state (1 for D*tt), and 5 is the strong 
phase difference between CFD and DCSD. Since the predicted value of R is small, ~ 0.02 [4], 
we neglect terms of O (R 2 ) (and hence take C = 1). The strong phase 5 for D*it is predicted 
to be small [3, 5]. Since R is expected to be suppressed, the amount of CP violation in D*tt 
decays, which is proportional to R, is expected to be small and a large data sample is needed 
in order to obtain sufficient sensitivity. We employ a partial reconstruction technique [6] 
for the D*tt analysis, wherein the signal is distinguished from background on the basis of 
kinematics of the 'fast' pion (iif) from the decay B — > D*7Tf, and the 'slow' pion ('7T S ') from 
the subsequent decay of D* — > Dn s ; the D meson is not reconstructed at all. 
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Previous analyses have been reported by Belle [7, 8] as well as by Babar [9]. This study 
uses a data sample of 605 ftT 1 containing 657 million BB events, which is about two times 
the size of the dataset used in the previous Belle analysis [8] and supersedes the previous 
Belle result. 

BELLE DETECTOR 

The data were collected with the Belle detector [10] at the KEKB asymmetric energy 
electron-positron (e~e + ) collider [11] operating at the T(4S) resonance. The Belle detec- 
tor is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel-like arrangement of time-of- flight scintillation counters (TOF), 
and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoidal coil that provides a 1.5 T magnetic field. An iron flux- return lo- 
cated outside of the coil is instrumented to detect K® mesons and to identify muons (KLM). 
A sample containing 152 million BB pairs was collected with a 2.0 cm radius beampipe 
and a 3-layer silicon vertex detector (SVD1), while a sample of 505 million BB pairs was 
collected with a 1.5 cm radius beampipe, a 4-layer silicon vertex detector (SVD2), and a 
small-cell inner drift chamber [12]. 

ANALYSIS PROCEDURE 

Partial Reconstruction of B — > D^tt^ decays 

To reconstruct the CP-side tags, we use: B —> P* =l= 7r =F ; D*^ — > D°7r ± . Candidate events 
are selected by requiring the presence of oppositely charged '/' an d V candidates. We 
estimate the D* frame using energy-momentum conservation: 

E D * = Ego — Ef, 

Pf+PD*=PB- (3) 

Here, E and p stand for energy and momentum respectively. E-^ is half the total centre- 
of-mass energy {Eqm) of the incoming e + e~ beams and ps = a/ E^ M /A — m' 2 B0 , where m B o 
is the nominal B° mass [14]. Using Ed*, the momentum of D* can be obtained in the e + e~ 
centre-of-mass (cms) as: Pd* = a/ Ej^ — m 2 D , . We construct a partially reconstructed D* + 
frame, using pr>* and Ed* and taking the direction of pn* opposite to pf. 

We define a variable, ps, which is strongly correlated with the fast pion momentum, pf. 
Ps is defined as: \\pj\ — \pd*\\ an d from Eq. (3), it follows: \ps\ < \pb\ (~ 0.3 GeV/c). 
We boost the charged slow pion into the partially reconstructed D* frame. In the true 
D* frame, the slow pion is mono-energetic. However, in the partially reconstructed D* + 
frame, the slow pion momentum will have a limited spread. We study the parallel and the 
transverse components of the momentum of the slow pion, ir s along the direction opposite 
to /, which are denoted as p\\ and p±, respectively. In the true D* frame, the p\\ variable 
has a distribution proportional to cos 9 2 for signal events, as the B decay is a pseudoscalar 
to pseudo-scalar vector transition. 
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CP side selection 

Fast pion (/) candidates are required to have a radial (longitudinal) impact parameter 
dr < 0.1cm (\dz\ < 2.0cm) and to have associated hits in the SVD. We reject leptons and 
kaons based on information from the CDC, TOF and ACC from the fast pion candidate list. 
A requirement is made on the fast pion cms momentum, 1.93GeV/c < pf < 2.50GeV/c. 
Slow pion (s) candidates are required to have cms momentum in the range 0.05GeV/c < 
p s < 0.30GeV/c. Since slow pions are not used for vertexing, no particle identification 
requirement is applied, we impose only a loose requirement that they originate from the IP. 
We select candidates that satisfy +0.00GeV/c < p± < +0.06GeV/c, -O.lOGeV/c < p\\ < 
0.07 GeV/c and -0.60 GeV/c < p s < 0.50 GeV/c . 

Vertexing and Flavour Tagging 

The determination of the flavour of the B meson opposite to the CP-side B is essential for 
the At measurement. In order to tag the flavour of the associated B meson, we require the 
presence of a high- momentum lepton (/) in the event. This helps reduce background from 
continuum e + e~ — > qq (q = u, d, s, c) processes. Tagging lepton candidates are required to 
be positively identified either as electrons, on the basis of information from the CDC, ECL 
and ACC, or as muons, on the basis of information from the CDC and the KLM. They are 
required to have momenta in the range 1.1 GeV/c < pi < 2.3 GeV/c, and to have an angle 
with the fast pion candidate that satisfies —0.75 < cos6 nf i in the cms. The lower bound on 
the momentum and the requirement on the angle also reduce, to a negligible level (0.7%), 
the contribution of leptons produced from semi-leptonic decays of the unreconstructed D 
mesons in the B° — > D*^^ decay chain. 

Identical vertexing requirements to those for fast pion candidates are made in order to 
obtain an accurate z tag position. To further suppress the small remaining continuum back- 
ground, we impose a loose requirement on the ratio of the second to zeroth Fox- Wolfram [13] 
moments, i?2 < 0.6. 

At the KEKB asymmetric-energy e + e~ (3.5 GeV on 8 GeV) collider, operating at the 
T(4S') resonance (a/s = 10.58 GeV), the T(4S') is produced with a Lorentz boost of (3^ = 
0.425, almost along the electron beamline (z) at KEKB. In the cms, B° and B° mesons 
are approximately at rest. Hence the proper time-difference (At) between the zqp and z ta . g 
vertices is obtained from fast pion on the CP-side and the tagging lepton. The variable At 
is defined as: 

At « (z CP - Ztoz)lfac. (4) 

The CP-side (zcp) vertex is obtained from the fast pion on the CP-side and the run- 
dependent interaction point profile (IP). The tag-side (ztag) vertex is obtained from tagging 
lepton and the run-dependent IP. 

Yield Fit 

We use the three kinematic variables, p§, p\\ and p± to distinguish between signal and 
background on the CP side. Background events are separated into three categories: P* T p ± , 
which is kinematically similar to the signal; correlated background, in which the slow pion 
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originates from the decay of a D* that originates from the decay of the same B as the 
fast pion candidate (e.g., D**n); and uncorrelated background, which includes everything 
else (e.g., continuum processes, Dtt). The kinematic distributions of the signal and back- 
ground categories are determined from a large MC sample, corresponding to three times the 
integrated luminosity of our data sample. 

We select candidates that satisfy -O.lOGeV/c < p\\ < +0.07GeV/c and -0.60GeV/c < 
Ps < +0.50GeV/c. In the cases where more than one candidate satisfies these criteria, 
we select the one with the largest value of 5 7rf7ts , where S nf7TB is the angle between the fast 
pion direction and the slow pion direction in the cms frame. The signal region is defined as: 
-0.40 GeV/c < p s < +0.40 GeV/c and two regions in pf -0.05 GeV/c < p\\ < -0.01 GeV/c 
and +0.01 GeV/c < p\\ < +0.04GeV/c. 

Event-by-event signal and background fractions are determined from binned maximum 
likelihood fits to the two-dimensional kinematic distributions (p$ and p»). The results of 
these fits, projected onto each of the two variables, are shown in Fig. 2, and summarized in 
Table I. We obtain a signal purity of 59.0 ± 0.4%, where purity is defined as ratio of the 
signal and total yield. 




FIG. 2: Results of the yield fits to D*-k candidates projected onto the p\\ (left) and ps (right) axes in 
the signal region of the kinematic variables. The contributions are: D*ir (yellow), D*p (magenta), 
correlated background (blue) and uncorrelated background (red). 



Fit Procedure to obtain CP violation parameters 

We perform a simultaneous unbinned fit to the same-flavour (SF) events, in which the fast 
pion and the tagging lepton have the same charge, and opposite-flavour (OF) events, in which 
the fast pion (/) and the tagging lepton (/) have the opposite charge, to measure the CP 
violation parameters in the D*tt sample. We minimize the quantity — 21n£ = — 2^1n£j, 
where 

C*i f D*ttPd*tt + f D* pPD* p + /unco -Punco + ./corr-fcorr- (5) 
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TABLE I: Summary of the yields in the signal region 



Candidates 

D*tt 50196 ± 286 

D*p 10232 ± 150 

Correlated background 10425 ± 135 
Uncorrelated background 14193 ± 128 



Here, / stands for the event-by-event signal and background fractions and are obtained from 
the fits to the kinematic variables and P stands for the probability density functions (PDF) 
for signal and backgrounds, which contain a physics PDF and experimental effects. For D*ir 
and D*p, the PDF is given by Eq. (1), where for D* p the terms S ± are effective parameters 
averaged over the helicity states [15] and are constrained to be zero. The PDF for correlated 
background contains a term for neutral B decays (given by Eq. (1) with S ± = 0), and a term 
for charged B decays (for which the PDF is 2 ^ e~\ At ^ T B+ , where rg+ is the lifetime of the 

charged B meson). The PDF for uncorrelated background also contains neutral and charged 
B components, with the remainder from continuum e + e~ — > qq (q = u, d, s, c) processes. 
The continuum PDF is modelled with two components: one with negligible lifetime, and the 
other with a finite lifetime. 

The parameters in P unco and P C orr are obtained from separate simultaneous fits to OF and 
SF candidates in the respective sideband regions. In these fits, the CP violation parameters 
are fixed to since there is no CP in background. The fit is further simplified by fixing the 
biases in Az to zero (discussed later in detail). Monte-Carlo simulation studies demonstrate 
that floating or fixing these biases to does not affect the background parameters. 

To measure the uncorrelated background shape, we use events in a sideband region, 
-O.lOGeV/c < p|| < 0.07GeV/c, 0.01 GeV/c < p\\ < 0.04GeV/c, -0.60GeV/c < p s < 
0.50 GeV/c and 0.08GeV/c < p± < O.lOGeV/c, which is populated mostly by uncorre- 
lated background. To determine the correlated background parameters, we use events in 
a sideband region, -O.lOGeV/c < p\\ < 0.07GeV/c, -0.60GeV/c < p s < O.OOGeV/c and 
O.OOGeV/c < p± < 0.05GeV/c. This sideband region is dominated by correlated and un- 
correlated backgrounds and has very small amount of D*ir signal and D*p background. The 
uncorrelated background parameters are fixed to the values obtained in the previous fit. 

The PDF for the signal and background in Eq. (5) must be convolved with corresponding 
Az resolution functions related to kinematic smearing (JZ^), detector resolution (JZdet), 
and asymmetry in Az due to non-primary tracks (lZ np ). The resolution function related to 
kinematic smearing is due to the fact that we use the approximation of Eq. (4). The detector 
resolution function parameters are obtained using J/ip — > /i 4 pT candidates. Since both the 
fast pion and the tagging lepton originate directly from B meson decays for correctly tagged 
signal events, we do not include any additional smearing due to non-primary tracks in such 
events. However, for incorrectly tagged events, almost exclusively originating from secondary 
leptons or pions, the PDF is convolved with an additional resolution component whose 
parameters are determined from MC simulations. The detector resolution and smearing due 
to asymmetry in Az due to non-primary tracks are described in detail elsewhere [8]. 
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Mistagging is taken into account using 

P(/ t =F g; vr±) = (i _ w T )P(B°/B° -> D*^) + w±P(B°/B° -> D* T 7r ± ) (6) 

where 717 is fast pion in CP-side, Z is tag-side lepton, u> + and w~ are the wrong-tag fractions, 
defined as the probabilities to incorrectly measure the flavour of tagging B° and B° mesons 
respectively and are determined from the data as free parameters in the fit for S^. 

The time difference At is related to the measured quantity Az as described in Eq. (4), 
with an additional term due to possible offsets in the mean value of Az, 

At — >At + e At ^ (Az + e Az ) I he. (7) 

It is essential to allow non-zero values of e At since a small bias can mimic the effect of CP 
violation: 

cos(AmAt) — > cos(AmAt) — Ame^t sin(AmAt) (8) 

A bias as small as e Az ~ 1 /im can lead to sine-like terms as large as 0.01, comparable to 
the expected size of the CP violation effect. Because both vertex positions are obtained 
from single tracks, the partial reconstruction analysis is more susceptible than other Belle 
CP analyses to such biases. We allow separate offsets for each combination of h and / 
charges. In order to correct for a known bias due to the relative misalignment of the SVD 
and CDC in SVD1 data, a small correction is applied to each measured vertex position. 
This correction is dependent on the track charge, momentum and polar angle, measured in 
the laboratory frame and is obtained by comparing the vertex positions calculated with the 
alignment constants used in the data, to those obtained with an improved set of alignment 
constants [16], which removes the observed bias. Since the alignment in SVD2 data was 
found to be comparable to that of the corrected SVD1 data, no additional correction was 
applied to SVD2 data. 



Fit Result 

In order to test our fit procedure, we first constrain S + and S~ to be zero and perform 
a fit in which r^o and Am (as well as two wrong tag fractions and eight offsets) are free 
parameters. We obtain r^o = 1.538 ± 0.008 ps and Am = 0.482 ± 0.004 ps" 1 , where the 
errors are statistical only. These values are compatible with their world average values [14] . 
Reasonable agreement with the input values is also obtained in MC. Furthermore, fits to 
the MC with S ± floated give results consistent with zero, as expected. 

To extract the CP violation parameters we fix r^o and Am at their world average values, 
and fit with S + , S~, two wrong tag fractions, and eight offsets as free parameters. We 
obtain S + = 0.057 ± 0.019 and S~ = 0.038 ± 0.020 where the errors are statistical only. 
The wrong tag fractions are W- = (6.8 ± 0.3)% and u>+ = (6.6 ± 0.3)%. All floating offsets 
are consistent with zero except for one of the OF combinations (h = 7r~, I = £ + ) in the 
SVD1 sample. The results are shown in Fig. 3. To further illustrate the CP violation effect, 
we define asymmetries in the same flavour events (v4. SF ) and in the opposite flavour events 
(A OF ), as 

(SF N n -i-(Az) - N 7V +i+(Az) 



A" 



.of N w+l -(Az) - N n - l+ (Az) 

' N n+l -(Az) + N 7T - l+ (Az)' {) 
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where the N values indicate the number of events for each combination of h and / charge. 
These are shown in Fig. 4. 




FIG. 3: Az distributions for 4 flavour-charge combinations: ir~l~ (top left) , ir~l + (top right), 
7r + /~ (bottom left), ir + l + (bottom right). The fit result is superimposed on the data (blue line). 
The signal and background components are shown as the red and dotted black curves, respectively. 



SF Asymmetry 

U.3 




"°-0\2-0.1 5-0.1-0.05 -0 0.05 0.1 0.15 0. 

Az(cm) 



OF Asymmetry 

U.3 




$.2-0. 15-0.1-0.05 -0 0.05 0.1 0.15 0. 

Az(cm) 



FIG. 4: Results of the fit to obtain S + and S , shown as asymmetries in the same flavour events 
(left) and opposite flavour events (right). The fit result (red curve) is superimposed on the data. 
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Systematic Error 

This analysis is very sensitive to the vertexing bias. Hence, we have use Az offsets in 
the fits to take care of this bias. In order to estimate the error due to these offsets, we use 
the difference of the mean of S ± obtained using an ensemble of 300 generated D*ir signal 
samples with CP (S ± = —0.04) and the generated CP value. 

Other sources of systematic error are the parameters of resolution functions, TZ^, TZ^ et 
and TZ np , the parameters of uncorrelated and correlated background and physics parameters, 
Am, tbo, t b +, Sp* and S^ rr that are fixed in the CP fit, where S^ rr are the CP viola- 
tion parameters for the correlated background component (S^ TI = ±0.05 in the CP fit). 
Additional systematic errors can result from varying the number of bins for the kinematic 
variables, ps and pii in the yield fit. 

We use a triple Gaussian as the detector resolution (i?det) function model. We consider 
the systematic uncertainty due to lack of knowledge of the exact functional form of the 
resolution model. Hence, we change the resolution models and obtain shifts as large as 0.008 
for S + . This is conservatively assigned as the systematic error due to lack of knowledge of 
the resolution model. 

We also performed a linearity test to check for possible fit bias by generating a number 
of large samples of signal MC simulations with different input values of S + and S~. All 
results are consistent with the input values, without evidence of any bias. In addition, we 
checked the pull for S + and S~ using two types of ensembles, one set generated with no 
CP (S^ = 0) and the other with CP (S ± = —0.04) and obtained mean (m) and rms (a) of 
the pull distributions fitted to a single Gaussian. For the no CP case, m s + = +0.10 ± 0.06, 
a s+ = +0.94 ± 0.05; m s - = -0.10 ± 0.06, a s - = +0.96 ± 0.05 and for the case with CP, 
m s+ = +0.14 ± 0.07, a s+ = +1.10 ± 0.06; m s - = -0.29 ± 0.06, a s - = +0.98 ± 0.05. Both 
cases yield m and a for the pull distributions close to and 1, respectively. This shows that 
our fit routine does not have any significant bias. 

The systematic errors are summarized in Table II. The total systematic error is obtained 
by adding the above terms in quadrature. 



Systematic error source 


S+ 


5" 


Az offset 


0.002 


0.003 


TZk parameters 


0.002 


0.003 


TZdet parameters 


0.002 


0.002 


TZ np parameters 


0.008 


0.007 


Background parameters 


0.002 


0.003 


Physics parameters 


0.004 


0.004 


Yield fit 


0.003 


0.003 


Resolution model 


0.006 


0.002 


Az floated in background PDF 0.000 


0.000 


Total systematic error 


0.012 0.010 



TABLE II: Summary of possible sources of systematic error 
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The results using the partial reconstruction method are 

S+ = +0.057 ±0.019 ±0.012, 

S- = +0.038 ±0.020 ±0.010, (10) 
where the first error is statistical and the second error is systematic. 

SUMMARY 

We have measured CP violation parameters that depend on 3 using the time- dependent 
decay rates of the decay B° — > _D* =F 7r =t: using a data sample containing 657 million BB 
events. We obtain the CP violation parameters expressed in terms of S + and S~, which are 
related to the CKM angles 4>\ and 03, the ratio of suppressed to favoured amplitudes, and 
the strong phase difference between them, as S ± = —Rd*-k sin(2</>i + <f) 3 ± 5d*-k)/ (1 + Rd*-k) 
for D*n as 

S + = +0.057 ±0.019 ±0.012, 

S- = +0.038 ±0.020 ±0.010, (11) 
where the first errors are statistical and the second errors are systematic. 
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